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\ Abstract: In this paper, an opto-thermal analysis of a moderately heated 

■ Ughtweighted solar telescope mirror is carried out using 3D finite element 

analysis (FEA). A physically realistic heat transfer model is developed to 
account for the radiative heating and energy exchange of the mirror with 
surroundings. The numerical simulations show the non-uniform temperature 
distribution and associated thermo-elastic distortions of the mirror blank 
clearly mimicking the underlying discrete geometry of the lightweighted 
substrate. The computed mechanical deformation data is analyzed with 
surface polynomials and the optical quality of the mirror is evaluated with 
(-H . the help of a ray-tracing software. The thermal print- through distortions 

pi are further shown to contribute to optical figure changes and mid-spatial 

^ i frequency errors of the mirror surface. A comparative study presented 

^ ' for three commonly used substrate materials, namely, Zerodur, Pyrex and 

i Silicon Carbide (SiC) is relevant to vast area of large optics requirements in 

' ground and space applications. 
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1. Introduction 

Large telescopes provide higher spatial resolution and improved sensitivity to probe exotic and 
unknown physical processes taking place in the Universe. Since last four hundred years, the 
telescope aperture size has gradually increased, doubling almost in every 50 years period IT]. 
For night time astronomy, this impressive growth will continue in the new millennia as the work 
on few notable 20-40 m class giant telescope projects is already in progress |l2]. The evolution 
of solar telescope, on the other hand, has followed a different trajectory ||3l. The immense dif- 
ficulty in managing the thermal loads and internal seeing caused by radiative flux has been the 
major hurdles in building large solar telescopes. The technology has now advanced to a stage 
where building large aperture solar telescopes is well within the reach. Efforts from various 
groups are now underway to build next generation 2-4 m class solar telescopes. Some notable 
solar facilities that are currently in different phases of development include: a 4m Advanced 
Technology Solar Telescope (ATST) to be build by the National Solar Observatory, USA I?]. 
The 4m European Solar Telescope (EST), under considerations from European Association for 
Solar Telescopes fS). A 2m class National Large Solar Telescope (NLST) proposed by the In- 
dian Institute of Astrophysics, India ||6l. These telescopes will be equipped with adaptive optics 



and other state-of-the-art scientific instruments to facilitate diffraction-limited observations of 
the small scale features (~50 km) of the solar atmosphere. Meanwhile, the 1.6 m New Solar 
Telescope (NST) at Big Bear Solar Observatory, USA and the 1.5 m Gregor telescope built 
by a consortia led Kiepenheuer-Institute, Germany, have started providing high quality ground 
based observations of the sun QISl. 

Evidently, one of the most challenging tasks is to carefully handle the excessive heat gener- 
ated by large radiation flux incident on the primary mirror Equally important is the choice of 
high quality, thermally stable and rigid material capable of maintaining the optical figure over a 
specified range of observing conditions. The thermal response of the solar mirror is influenced 
by radiative heating of the mirror caused by direct light absorption and the temperature im- 
balance with the surroundings. Mirror heating adversely affects the final image quality of the 
telescope in two ways, as is explained below. 

First and foremost, it is hard for a massive mirror having large heat capacity and high ther- 
mal resistance to attain a temperature balance with the surrounding air Under such conditions, 
the temperature difference between the mirror surface and the ambient air leads to refractive 
index fluctuations along the beam path. This 'mirror seeing' produces wavefront aberrations 
responsible for the image blurring. To preserve the image quality of the telescope, the effect of 
mirror seeing has to be minimized. This means the temperature difference between the mirror 
surface and the ambient air should not exceed beyond ±2° C ||9|. More crucially, the tem- 
perature uniformity across the mirror surface has to be maintained within ±0.5° C. Several 
approaches which include, air conditioning ifTOl . fanning the optical surface ifTTI . ventilating 
mirror interior, resistive heating [112JI131 have been proposed to regulate the mirror temperature 
for ensuring high quality astronomical observations. Laminar air flow across the mirror ho- 
mogenizes the temperature by disintegrating the 'thermal plumes' closer to the surface. These 
methods also speed up the thermalisation process by improving the convective heat transfer rate 
between the substrate and the surrounding air 

Secondly, a poorly conducting mirror is differentially heated to create temperature inhomo- 
geneities within the substrate volume. Thermally induced material expansion or contraction 
can significantly alters the optical quality of the telescope mirror This problem can be elimi- 
nated by choosing mirror substrates made of ultra low expansion (ULE) material. In fact, the 
demand for high quality astronomical optics has led to the development of many attractive low 
expansion materials such as Zerodur (Schott), fused quartz, Cervit, ULE (Coming), Clearce- 
ram (Ohara) and Astro-Sitall etc. These ULE materials usually have high heat capacity and low 
thermal conductivity, implying they can retain a significant amount of heat, leaving the mirror 
prone to undesirable seeing effects and surface buckling under high thermal loads. 

The severity of mirror heating problem scales adversely with the aperture size. To reduce 
the overall weight and large thermal inertia, the most preferred choice is to use lightweighted 
mirror structures that are easy to mount and do not deform significantly under gravity. Besides, 
the reduced mass also facilitates efficient and faster cooling of the mirror in harsh thermal 
conditions. For space applications, reducing the weight of the mirror has greater premium on 
the overall cost reduction of the payload. 

The lightweight geometry is created by removing pockets of material from the blank with- 
out significantly compromising the rigidity and stiffness of the mirror |[T4l . The extreme 
lightweighting is usually achieved by mechanical milling and acid etching. The cool air, cir- 
culated through the rear pockets can remove the excessive heat from the mirror very efficiently. 
The pocketed cell geometry introduces structural inhomogeneities within the mirror blank. The 
top faceplate (reflecting surface) of the mirror gets the mechanical support from the grid of 
discrete ribs. The effect of underlying geometry of the cellular structure usually manifests in 
the form of rib print-through patterns a) during the mechanical polishing of the faceplate and b) 



periodic temperature residuals on the mirror faceplate operating in extreme thermal conditions. 
While polishing, the lapping tool experiences greater reaction force from the rib locations to 
wear out more materials thus creating periodic thickness variations on the faceplate. The wave- 
front errors related to polishing print-through are well studied and there are ways to mitigate 
them IBHnl. 

In this paper we have examined the impact of relatively less explored but nonetheless im- 
portant 'thermal print-through' issue associated with the lightweighted structure operating in 
thermally dynamic environment. The thermo-optic analysis, as it is increasingly recognized, is 
a vital component of the overall system design and necessary to predict and verify the satisfac- 
tory performance of the mirror [118JI191 . Here, we have used FEA to solve the time dependent 
3-D heat transfer equation to predict the mirror deformation arising from the thermally induced 
stress and temperature inhomogeneities of a lightweighted structure. For realistic numerical 
simulations, the location dependent solar flux and ambient heating model is incorporated into 
FEA analysis. Thermal response of three types of mirror materials is studied at varying am- 
bient temperature conditions that may exist at various observatory locations. Temperature in- 
duced surface deformations were analyzed using Zemike polynomials. For optical analysis, the 
Zemike data was imported into Zemax -a ray tracing software. A comparative study is carried 
out for three well accepted substrate materials in astronomical optics, namely, Zerodur, Pyrex 
and SiC. 

2. Mathematical formulation 

In order to study the primary mirror heating and the associated thermal stress/strain fields af- 
fecting the surface figure, each physical process should be identified. In the following we give a 
brief mathematical formulation of the heat transfer problem along with appropriate prescription 
for boundary and initial conditions. 

2.1. Heat transfer 

In a simplified form, the telescope mirror can be considered as a 3D circular disk whose front 
face is exposed to intense solar irradiance /(W/m^). The thin metallic coating on the faceplate 
absorbs about 10%-15% of radiant flux which is converted into heat. The temperature distribu- 
tion inside the mirror changes as the heat conducts from front surface to the rear side. The 
mirror also exchanges heat (via convection and radiation) with ambient air, which is at a tem- 
perature Ta. A heat transfer equation can be derived based on the principle of conservation of 
energy, solution of which gives the temperature distribution r(x,y,z,f). The 2nd order partial 
differential equation relevant to the current studies for the conductive heat transfer in a solid 
can be expressed as ll20l 

where T is the temperature, K = k/ (pCp) is the thermal diffusivity of the system, Cp is the heat 
capacity, k is the thermal conductivity, p is the density and q = Q/{pCp), where Q{x,y,z',t) 
represents the source term for heat. In the present case there is no internal heat source, i.e. 
Q{x,y,z',t) = everywhere and at all instants of time. The heat flux enter the system through 
the surfaces, a fact which must be accommodated in the boundary condition. 
For solving the temperature T{x,y,z,t) we must assign; 



- the initial condition 

T{x,y,z;t = 0) = Ti{x,y,x) and 



(2) 



- the boundary condition that the net incident flux, which is balance between the incident 
and emitted fluxes, is transmitted into the substrate and accounts for the temperature 
gradients on the domain surface Q.. This condition is described by, 

k (Vr • n) In = [^0 -hiT~ T^) - e^a^ {T^ - T^)] \a (3) 

where n is surface normal vector at any point on CI. In Eq. (3) the l.h.s. denotes the heat inflow 
due to temperature gradient on the surface, qo is incident flux due to solar radiation, h {T — 7a) 
describes the convective heat flow to surroundings, where 7a (f) is ambient temperature and h 
is the heat transfer coefficient. The term EsCTr (r"* — T^) describes the radiative exchange with 
surroundings, where Es is the emissivity of the surface and (7r is Stefan-Boltzmann constant. 

2.2. Model prescription for initial temperature, heat flux and ambient temperature 

2.2.1. Initial temperature Ti (x, y, z) 

Due to poor thermal conductivity and asymmetric heating, the initial temperature distribution 
rj(x,y,z) inside the mirror cannot be assumed to be uniform, i.e., independent of {x,y,z). In 
order to establish a definitive temperature field at f = 0, simulations were carried out by sub- 
jecting the mirror to several day-time heating and night-time cooling cycles ||2T1 . The resulting 
temperature distribution was used as initial 7/ for all subsequent studies. 

2.2.2. Surface boundary conditions for inward heat flux qQ 

The boundary prescription for inward heat flux in Eq. (3) depends on the solar irradiance 
7(W/m2) that varies between the sun-rise and the sun-set as lISTl 

qo = yl and (4) 
7 = 7o [sin((p)sin(5) +cos(9)cos(5)cos(77)] (5) 

where y is absorption coefficient (fraction of radiant energy that is absorbed and converted into 
heat) of the metallic film on the reflecting surface, 7o(W/m^) is irradiance amplitude when sun is 
at zenith, 5 (rad) is declination angle of the sun, (p (rad) is latitude of the place and H = {n/l2)t 
is solar hour angle and t is the local time. 
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Fig. 1. Typical variations in ambient temperature (red curve, scale: vertical-left) and solar 
flux reaching the mirror surface, (blue curve, scale: vertical-right). The peaks of solar flux 
and ambient temperature are shifted by 2 hrs due to thermal time lag. 



2.2.3. A simple model for ambient temperature Ta 



In Eq. (3), the heat loss or gain by free convection and radiation is driven by ambient 
temperature T^. For convenience, we can split the ambient temperature into two parts, i.e., 
7a (0 = 7b(0 + 7n(0' where TY){t) and Tfi{t) describe the day-time ambient heating and night- 
time cooling, respectively. For cloud-free sky conditions, these two terms can be approximated 
using a simple model given by Gottsche and Olesen 



7b(0 = 7;+7bcos 



0) 



for t < ts (6) 



Tn{t) = 7;. + rrf+{rocos [^(r-T,,,)] -7;,}exp[-(f-f,)/'f] for f>fv (7) 

We use Eq. (6) and Eq. (7) as input to our FEA model to simulate a range of ambient 
temperature conditions typically existing at observatory locations. One such realization for 
10°C < 7a(f) < 20°C and solar flux /(W/m') is shown in Fig. 1. The assumed values for var- 
ious parameters in Eq. (6) and Eq. (7) used to obtain a smoothly varying temperature curve in 
Fig. 1 are EH: To = 6.8°C; T, = 13.4°C; T,, = -3.5°C; CO = 12hrs (the half period oscillation); 
T„, = 14 : 00 hrs (the time when temperature maxima is reached); f j = 18 : 00 hrs (the time when 
temperature attenuation begins); K" = 3.5 hrs (the attenuation constant). The ambient tempera- 
ture 7a shown in Fig. 1 continues to rise well past the solar noon. This is due to thermal time 
lag in terrestrial environment that exists between the peak of solar flux and the peak of ambient 
temperature. 

2.3. Structural deformation 

The heat conducted to different parts of the mirror builds up thermal strain which results in 
structural deformation. For isotropic linear elastic solid, the thermal strain depends on coef- 
ficient of thermal expansion a(K^'), the instantaneous temperature T{t) and the stress-free 
reference temperature T^ef as: 

eth = a(r-7;-ef) (8) 

If M, V, and w are the thermally induced deformation components in x, y and z direction, then 
the thermal strain can be completely specified in terms of m, v, and w and their derivatives. For 
small displacements, the constitutive equations relating the components of normal and shear 
strain to the deformation derivatives can be expressed as 



(9) 



The displacements at the bottom boundary of the mirror were kept fully constrained, i.e., u 
v = w = 0, to avoid the rigid body motion. 
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3. The mirror geometry 

Among several possibilities of lightweighted geometry, the open-back mirror configuration is 
highly suitable for solar telescopes ||241 . It consists of a thin faceplate fused to an array of ribs 
formed by scooped-out material from the back. The pocket geometry can be circular, triangular. 
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Fig. 2. Mirror geometry: (a) The side view and (b) the CAD model of the mirror with open 
back lightweighted honeycomb structure and 12 mm thick faceplate on top. Mirror has 6- 
fold symmetiy about the z-axis. The faceplate is shown placed only on one of the 6 sectors 
of the mirror. 



hexagonal, square or any other combination of these shapes. The resulting structure has side 
walls/ribs of certain thickness that supports the thin reflecting faceplate of the mirror from 
behind. Besides the mechanical support, the rib walls also provide the thermal flow path for 
heat to conduct away from the faceplate. The presence of discrete grid makes the heat removal 
rate to vary across the faceplate thus creating temperature impressions resembling the cellular 
geometry of the lightweighted structure. 

Figure 2(a) shows the schematic of the parabolic mirror (conic constant =-1, focal length 
4m) geometry chosen for the present studies. A similar design is also under consideration for 
the proposed NLST by India. The diameter of the mirror is 2m and the diameter of the inner 
hole is 0.42m. The honeycomb bottom comprising about 708 hexagonal cells is enclosed by 
20 mm thick outer rim. The side length L of the hex-cell (inset image in Fig. 2(b)) is 34.6 mm 
and the wall/rib thickness f is 8 mm. The diameter D of the inscribed circle is 60 mm while 
the wall/rib depth varies from 148 mm for the outermost cells to 85 mm for the one at the 
center The faceplate thickness is 12 mm. Compared to a solid geometry, the overall mass of 
the current lightweighted mirror is reduced by ~ 68%. The main reason for choosing these 
hex-cell dimensions is because their structural stability has already been verified for the Gregor 
1.5 m mirror telescope ll25l . The mirror has a 6-fold symmetry about the optical axis as seen in 
the CAD design shown in Fig. 2(b). The symmetry property is utilized to reduce the model size 
and computational requirements by selecting only one section of the mirror in the FEA studies. 
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Mirror radius (m) 

Fig. 3. The impact of daytime heating on the Zerodur mirror with ambient temperature 
varying between 20-30° C and h = 5 W/m-. (a) Temperature print-through pattern on the 
mirror faceplate subjected to da5'time radiative heating, (b) Temperature variations along 
the radial line AA' of the mirror faceplate at different times. 

4. The FEA results 

The governing equations and associated boundary conditions for the heat transfer and the struc- 
tural deformation problem are outlined in Section 2. The analytical solution for Eqs. (l)-(3) 
and Eqs. (8)-(9) exists only for simple and regular shapes. In most practical cases where ob- 
ject geometry is highly complex, boundary conditions are complicated, material properties are 
temperature dependent and various physical parameters are coupled, problem then has to be 
solved by numerical methods. We used COMSOL Multiphysics as our FEA tool to solve the 
time-dependent 3D heat transfer and structural problem. The segregated solver in COMSOL 
takes advantage of the direct coupling between heat transfer and structural model. That is, for 
each time step (At = 100 sec), the heat transfer module first solves for the temperature distribu- 
tion of the mirror over 24 hours. Each temperature field is then passed onto the the structural 
analysis module to perform the stress-strain analysis. No data transformation was necessary 
since a common mesh is applicable to both thermal as well as structural models in COMSOL. 



Table 1. Values of material properties used in FEA studies 



Property 




Material 




Units 




Zerodur 


SiC 


Pyrex 


Al 




CTE(a) 


-0.08 


3.7 


3.25 


25 


lO^^K"' 


Thermal conductivity (k) 


1.46 


118 


1.14 


237 


WK^'m^i 


Heat capacity (Cp) 


821 


672 


750 


910 


Jkg-'K-i 


Density (p) 


2.53 


3.20 


2.23 


2.70 


lO^-kgm"^ 


Emissivity (e,,) 


0.9 


0.9 


0.95 


0.05 




Young's modulus (E) 


90.3 


414 


64 


70 


10*^ Pa 


Poisson's ratio (v) 


0.24 


0.14 


0.2 


0.33 





A 3D CAD model of the mirror conforming to optical prescriptions is shown in Fig. 2. The 
reflecting surface of the faceplate was modelled using a 100 nm thick aluminium coating on the 
faceplate. The FEA model was run on the 1 /6th axisymmetric section of the mirror that was 
partitioned into 100644 tetrahedral mesh elements. The top surface of the mirror was continu- 
ously heated by time varying irradiated solar flux plotted in Fig. 1. The thermal and structural 
response was evaluated at three different ranges of ambient temperature 7a i.e., 0°C— 10°C, 



10°C— 20°C and 20°C— 30°C. The front, back and the outer rim of the mirror were subjected to 
convective cooling with the heat transfer coefficient h = 5 W/m^ representing natural cooling 
and h = 15 W/m^ representing a moderately forced cooling. In addition, surface-to-ambient 
radiation condition was specified for the faceplate. Also the symmetry boundary condition i.e., 
k (Vr • n) =0, was imposed on the side edges of the mirror In the FEA model, the mirror 
heating was 'on' (Eq. (4)) between sunrise (6:00 hrs) and sunset (18:00 hrs) and 'off' for rest 
of the time. The standard values of the material parameters used in FEA studies are listed in 
Table I ED. 

4.1. Thermal print-through and temperature nonuniformities 

Figure 3(a) shows the temperature distribution within the mirror volume at 14:00 hrs when 
ambient temperature also reaches its peak. Apart from strong thermal gradients along axial di- 
rection, the appearance of periodic temperature print-through pattern on the faceplate is clearly 
evident. The temperature excess Ar = T — T^ between the mirror and the ambient air is plot- 
ted in Fig. 3(b). The temperature nonuniformities on the mirror surface are quantified by two 
key variations. First, the cyclic undulations 5T (the temperature print-through effect) along the 
radial line AA' which is caused by differential heat flow from the faceplate to the underlying 
hex-cell structure. Second, the trailing temperature (edge-effect) which denotes the overall 
temperature difference between the mirror surface and the edges. This difference arises due to 
more efficient cooUng from the walls of inner hole and the outer periphery of the mirror Ac- 
cording to data shown for Zerodur in Fig. 3(b), the £,T could be as high as 2°C at noon. From 
the temperature distribution alone, we can loosely identify £, T being associated with low order 
surface aberrations while 5T contributes to mid-spatial errors as explained in the next section. 
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Fig. 4. Temperature difference between two reference points located at hexagonal cell- 
center (C) and cell-edge (E). 

The maximal departure of AT^^^, ST^^x and (^Tmax for an ambient temperature range 7a : 
20°C— 30°C for Zerodur, Pyrex and SiC material is listed in Table 2. From the Table 2, we may 
note that the amplitude of both AT" and 8 T depend crucially on the thermal conductivity of 
the material and the convective heat transfer coefficients used in simulations. A more efficient 
cooling mechanism is necessary, especially in Zerodur and Pyrex, to remove the excessive heat 
to minimize the mirror seeing effects. 

A diurnal variations of 5T computed between two reference points located at the rib wall 
and the center of hex-cell is shown in Fig. 4. Even though the surface temperature is always 
above the ambient, the impact of Ar and 5T are manifestly more pronounced during the noon 



Table 2. The temperature difference for different glass material computed for heat transfer 
coefficient hi =5 W/m^ and /12 = 15 W/m-. 





Ta range: 20°C-30°C 


Material 


Ar 


nax 


ST 


nax 


^ Tmax 




In 


hi 


In 


hi 


hi 


hi 


Zerodur 


5.38 


2.56 


0.36 


0.20 


2.0 


0.9 


Pyrex 


6.07 


2.82 


0.43 


0.21 


2.0 


1.0 


SiC 


2.56 


1.02 


0.04 


0.02 


0.5 


0.25 




0.6 0.8 
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Fig. 5. Thermally induced structural deformation: (a) for the Zerodur mirror and the cor- 
responding z-displacements w (nm) along the radial line AA' of the mirror faceplate for (a) 
Zerodur, (b) for Pyrex and (c) for SiC at different times of the day. The ambient temperature 
range was 20-30° C and h = 5 W/m^. 



hours. 

4.2. Thermally induced surface deformation 

The non-uniform temperature field such as shown in Fig. 3(a), causes the mirror substrate to 
deviate away from its nominal defined shape. Thermally induced shape change can be highly 
irregular in lightweighted structures. The surface distortions for each temperature input field 
were computed by solving the FEM structural model. The reference temperature Tj-gf for the 
undeformed mirror geometry is assumed to be 20 °C. The impact of lightweighted hex-cell ge- 
ometry is clearly noticeable in the simulated image showing the surface deformation in Fig. 5(a) 
for Zerodur mirror The corresponding axial displacements w along the radial line A4' of the 
mirror faceplate for different substrate materials are also plotted in Figs. 5(b)-5(d). The surface 



displacement w is positive in Pyrex and SiC substrates and gradually increases in radially out- 
ward direction. However, for Zerodur material, w can be both positive or negative because of 
zero crossing of the CTE around 25 °C. For a given temperature field, the surface undulations, 
represented by 5w, are not uniform over the entire mirror surface but exhibit same periodicity 
as the underlying hex-cell structure. The displacement close to the mirror edges is expressed 
by ^w. The gradient across the mirror surface is approximated by the slope w' ~ dw/dr. The 
numerical values computed for 5w,^w and w' during the peak mirror temperature at f = 14 ; 00 
are summarized in Table 3. 



Table 3. Parameters describing the mechanical deformation of the surface. 
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Fig. 6. (a) Surface sag for thermally deformed Pyrex mirror at f = 14 hrs and h = 5 W/m^. 
The variation of RMS surface sag with daytime (b) for Pyrex and SiC and (c) for Zerodur 
mirror. 



The FEM studies were also performed at other temperature range of interest, i.e., 0°C— 10°C 
and 10°C— 20°C. The surface deformation data was found almost similar to that given in Table 
3. 

5. Optical analysis 

5.1. Surface sag and Zernike polynomials 

The thermally induced deformation data obtained from FEA in the previous section has to be 
transformed into a useful form to evaluate the optical performance of the mirror A convenient 
way to model the deformation of an optical surface is to convert the displacement data to sur- 
face sag values. The finite element computed node displacement of an arbitrary point on the 
mirror surface can be expressed 2& d = {u^ + + w^)'/^, where u, v and w are components 
of displacement vector d along x, y and z direction, respectively. For small perturbations, the 
surface sag change As and the node displacement at a radial point r„ can be related by the 



following formula 



As = 



dw{ro) 



(10) 







cos 


— sin0 





" 




X 


y 




sin 9 


cos 9 










y 


z' 










1 







z 


As' 













1 




As 



The sag data for the axisymmetric mirror segment was imported into MATLAB and 
converted into full aperture using co-ordinate transformation of the form: 



(11) 



where 9 =N n/3 and = 2,3,4,5 and 6. Equation(ll) maps the surface sag data from the 
original mirror segment (unprimed coordinate) to other (primed coordinate) sectors. 

The MATLAB function TriScatteredlnterp which makes use of Delaunay triangula- 
tion algorithm, was used for interpolating irregularly spaced sag data over an uniformly sam- 
pled array of points. A typical example of surface sag for full aperture Pyrex mirror is shown in 
Fig. 6(a). The sag values vary between 2.8/xm— 6.13/xm and the root mean square (RMS) sag 
over the full aperture is 5.3/xm. The RMS wavefront error is twice the RMS sag. The diurnal 
variation of RMS sag for Pyrex, SiC and Zerodur, shown in Figs. 6(b) and 6(c), is consistent 
with variation of their CTE under specified temperature range. The optical aberrations of the 



0.02 
0.00 
-0.02 



, | f . ifc i> ,1,, 



«rj- 



I Pyrex | 



E 0.000- 
-0.008- 



20 40 60 80 100 120 140 160 180 200 20 40 60 80 100 120 140 160 180 200 

0.008 T 



I r ■ 



I I 



(a) 



0.00 |i. I ■ , 
-0.01 



20 40 60 80 100 120 140 160 180 200 
Zernike terms 



20 40 60 80 100 120 140 160 180 200 20 40 60 80 100 120 140 160 180 200 

0.01 - 



,L i ll , i |, J 

(b) I ' 



20 40 60 80 100 120 140 160 180 200 

Zernike terms 



Fig. 7. Zernike expansion coefficients corresponding to (a) defocus, primary, secondary and 
other higher order spherical aberrations which contribute to overall all figure change and 
(b) the remaining low-amplitude terms which contribute to mid-frequency surface errors. 



thermally deformed surface were further analyzed by fitting Zernike annular polynomials to 2D 
sag data. The main advantages of Zernike polynomials arise due to their orthogonality over the 
unit circle, rotational invariance, completeness property and direct relationship of expansion 
coefficients with known aberrations of the optical system 027112811 . Each Zernike term is a prod- 
uct of three components, namely, a normalization factor, a radial part and an azimuthal part of 
type cosm(j) or sinm0. The surface sag As{p,(p;e) for an annular pupil with obscuration ratio 
e, radial variable p and azimuth angle can be written in terms of Zernike annular polynomial 
as: 

A5(p,0;e) = ^ fl;ZXp,0;e) (12) 

7=1 

where j is a single index polynomial-ordering term and aj are expansion coefficients ll29l . 




Fig. 8. Residual map obtained by subtracting the reconstructed Zemike surface from the 
original sag data. 

The j-th order Zemike polynomials Zj{p,(j);e) in Eq. (12) is defines as ll28l 

{^2(« + l)/?;"(p;e)cosm^ m 7^0,;' even 
v /2(« + l) (p ; e ) sin m 7^ 0, ; odd (13) 
yA^Rf,{p;e) m = Q 

where n is radial degree and m is azimuthal frequency such that m < n and n — \m\ is even. 
The procedure to obtain expressions for radial component iJ™ (p ; e ) in Eq. ( 1 3) is discussed in 
Ref Gl. 



Table 4. Low-order radial Zemike aberration terms. 
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We used MATLAB based open-source toolkit developed at the University of Arizona, for 
the Zernike analysis of the sag data 1301 . The computed Zernike coefficients for three substrate 
materials are plotted in a bar chart shown in Fig. 7. Based on their relative amplitude, various 
Zernike terms are separated into two parts. The leading contribution comes from purely radial 
part (i.e., m = case in Eq. (13)) plotted in Fig. 7(a). These terms correspond to piston (j = 1), 
defocus (j ~ 4), primary spherical (j ~ 11), secondary spherical {j = 22,) and higher order 
(7 = 37,56,79, 106, 137, 172) spherical aberrations that are responsible for driving the overall 
figure change of the mirror. Among these, the defocus term (Fig. 7) dominates for Zerodur, 
while primary spherical aberration dominates for Pyrex and SiC. The first 4 radial polynomials 
and their associated meanings are given in Table 4. The amplitude of these terms diminishes 
progressively with increasing order j. The pure radial and few low-frequency azimuthal terms 
are also reported to have strong influence on ocular aberrations encountered in vision related 
problems ||3T1 . 

The remaining low-amplitude Zemike terms are plotted in the Fig. 7(b). The resulting am- 
plitude of these mixed (product of radial and angular) Zemike terms is about 2-order smaller 



Fig. 9. Optical model of a Gregorian-type reflecting telescope in Zemax. Ml: primary mir- 
ror; M2: secondary mirror; Fl: primary focus; F2: secondary focus or image plane. 



Table 5. Predominant low-amplitude Zemike terms consisting of azimuth and radial parts. 
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than the pure radial terms. The predominant low-ampUtude terms comprising radial and angular 
part are given in Table 5. Ordinarily the higher order terms are not identified with any specific 
optical aberration. The solitary presence of cos 60 and its harmonics can be associated with 
the shape of thermal print-through errors arising from the 6 -fold symmetry of the hexagonal 
pockets. It is to be noted that amplitude of these terms (see Fig. 7(b)), does not decline appre- 
ciably with increasing j. This implies that using more number of Zemike terms is unlikely to 
improve the fitting accuracy of the sag data further One reason for such behavior is the limited 
ability of the Zernike polynomials to capture high-frequency surface imperfections as oppose 
to their effectiveness in capturing low-frequency global changes in the surface shape. That is 
why in Zernike representation, the impact of high frequency thermal print-through on surface 
sag is likely to be swamped by low-order figure changes induced by defocus and all spherical 
terms. In order to effectively visualize small-scale surface elevations, we computed residual 
sag by extracting the Zemike reconstructed surface from the original sag data. The mid-spatial 
frequency features, typical of thermal-print through affects, are clearly seen in the residual map 
shown in Fig. 8. 




Fig. 10. The spot diagram showing the telescope performance for a on-axis 0° field point. 
The blue traces represent the ray distribution in the image plane due to temperature in- 
duced surface imperfections in Ml. The pattern has an envelop of six-fold symmetry due 
to hexagonal nature of the elements, as expected. The diffraction limited Airy disc is indi- 
cated by the black circle. The specified wavelength is 0.55 flm and the radius of the airy 
disc is close to 15 ^m. Top row (a)-(c), when all the Zemike terms were included. Bottom 
row (d)-(f), without piston, focus and spherical terms. 



5.2. Image quality 

The surface aberration data of thermally deformed mirror was analyzed using Zernike polyno- 
mials and the optical performance was evaluated in a commercial ray tracing software Zemax. 
A simple telescope design was created in Zemax as shown in Fig. 9. The design conforms to 
optical prescription of Ml given in Section 3. The aberration data can be incorporated into the 
optical model of the telescope by assigning a suitable surface type -the Zernike Standard Sag, 
to the Ml in Zemax. This surface type is commonly used for describing mechanical deforma- 
tions. Selecting Zemike Standard Sag surface allows user to define the standard surface (e.g. 
plane, spheres, conies etc) plus additional Zemike terms that can be imposed on the standard 
surface to represent surface aberrations [l32l . The aberration data for Ml along with number 
of Zernike terms, normalized radius of the aperture and a list of externally computed Zemike 
coefficients was directly imported into Zemax using Extra Data Editor The M2 was assumed 
to be free from any surface irregularity. 

In Zemax, there are many ways to assess the optical performance of the system. We used 
spot diagrams and modulation transfer function to evaluate the performance of the telescope 
system. The thermal distortions are assumed to become worse when the mirror temperature 
reaches its maximum at f = 14 : 00 hrs. For the worst case scenario, the effect of thermally in- 
duced aberrations on the telescope image quality can be readily seen in the spot diagram shown 
in Fig. 10. Except for Zerodur material, the RMS spot radius (Figs. 10(a)-10(c)) for Pyrex and 
SiC is significantly larger than the size of the airy disc, implying the diffraction limited perfor- 
mance of Ml made of high CTE materials, is severely compromised by thermal heating. The 
contribution of low amplitude Zemike terms (Fig. 7(b)) is estimated by setting the radial terms 
(Fig. 7(a)) to zero in Zemax. Clearly, the low amplitude terms add significant scatter to the en- 
ergy distribution as shown in Figs. 10(d)-10(f). The envelop of the spot scatter is also indicative 



of the periodicity and shape of the 'surface bumps' caused by uneven temperature distribution 
resulting from underlying hex-cell structure of the Ml. Another way to examine the optical 
performance is to compute the contrast variations of the imaging system at different spatial 
frequencies. The variation of contrast, i.e., MTF for Pyrex and SiC is plotted in Figs. 1 1(a) and 
11(b), respectively. Considering the severity of the thermal distortions, a steep fall in contrast 
at relatively low-frequencies is not completely unexpected. 

The preceding studies were also repeated for convective heat transfer coefficient h = 
15 W/m^. The corresponding RMS spot radii for the Zerodur, Pyrex and SiC were found to be 
3.3/im, 1305 /im and 1798/im. This is only a moderate improvement. This shows that forced 
circulation of ambient air alone in high CTE materials (Pyrex and SiC) is not enough to keep 
the thermal distortions within tolerable limits. For effective thermal control, the temperature of 
the circulated air has to regulated few degrees below the ambient temperature. A detailed CFD 
analysis is however needed to arrive at the optimal coolant air temperature that would minimize 
the mirror seeing and thermal distortions of the Ml. The numerical studies presented in this pa- 
per would complement those efforts. Finally, we draw the attention of the scientific community 
to consider the present scheme of opto-thermal analysis for optimization of the solar telescope 
design. 
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Fig. 11. The system contrast in terms of transverse modulation transfer function at various 
spatial frequencies for a 0° field angle (a) for Pyrex and (b) for the SiC mirror. For Zerodur 
(not shown), the ray aberrations are well within the diffraction limits and three MTF curves 
would overlap. 





6. Conclusions 

The temperature stabilization of a primary mirror system is one of the most challenging tasks in 
building a large aperture solar telescope. To achieve the desired optical performance of the tele- 
scope, the thermal response of the primary mirror needs to be accurately predicted and suitably 
controlled under the harsh observing conditions. For large aperture telescope, the lightweighted 
mirror offer significant advantages over traditionally used solid blanks. However, the pocketed 
cell structure created by the lightweighting process leads to geometry dependent material inho- 
mogeneities and non-uniform temperature patterns. Besides providing physical support to the 



faceplate, the rib walls of the lightweighted structure also serve as discrete channels for the heat 
flow, giving rise to periodic thermal print-through errors across the mirror surface. We carried 
of a detailed thermo-optic analysis of a 2m class solar telescope mirror A time dependent heat 
transfer and structural model of the lightweighted mirror was solved using FEA methods. The 
mirror faceplate was heated by a radiant energy from the sun. The distinct appearance of ther- 
mal print-through patterns markedly depends on the thickness of the faceplate and changing 
thermal environment of the surroundings. The thermal displacement FEA data was converted 
to surface sag values useful for evaluating the optical quality of the telescope. The Zernike anal- 
ysis of the sag data showed that while the major contribution to the optical figure error comes 
from primary, secondary and higher order spherical aberration terms, the thermal print-through 
arising from underlying hex-cell geometry mainly contributes to low amplitude, mid-spatial 
frequency surface errors. A high scatter in ray-traced RMS spot size and a steep fall in MTF 
signifies a severe degradation in optical image quality of high CTE (SiC and Pyrex) materials. 
While SiC and Pyrex materials show a greater propensity for thermal distortions. But that alone 
may not be a sufficient ground to completely rule out their usages at least in moderate sized 
mirrors. It is also true that SiC and pyrex may not be able to match the performance level of 
ULE materials, but in principle, the structural instability can be better ensured by designing an 
efficient temperature control system. A higher convective heat transfer coefficient showed only 
a marginal improvement in the image quality which also reinforces the need for a better cool- 
ing mechanism to remove the excessive heat from the Ml. The approach outlined in this paper 
would be useful to design an effective thermal control and temperature stabilization system for 
solar observations. 
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